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Abstract 

. We propose a simple one sided Monte-Carlo algorithm to distinguish, to any given degree of 

certainty, the symplectic group C„(q) = PSp2„((7) from the orthogonal group Bn{q) = f22n+i(g) 
where g > 3 is odd and n and q are given. The algorithm does not use an order oracle and works in 
(~| ■ polynomial, of nlogg, time. 

. This paper corrects an error in the previously published version of the algorithm [Q]. 
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^ . 1 Introduction 

, The aim of this paper is to correct an error in the previously published version Q of a Monte-Carlo 

' algorithm which distinguishes, to any given degree of certainty, the symplectic group C„(g) = PSp2„((j') 

, from the orthogonal group i?„(q) = 02„+i(g) where g > 3 is odd and n and q are given. Namely we 

prove the following theorem. 

I Theorem 1.1 There is a polynomial time one sided Monte- Carlo algorithm which, when given a black 

■ box group G isomorphic to a simple classical group of type or Cn defined over some field of given 

size q > 3 and odd characteristic, can decide whether G is isomorphic to Gn or not. If G is Gn then, 
with probability > 1/2, the algorithm produces a witness against isomorphism G ~ Bn{q). 



B 

> 

; — r ^ o 

I who corrected a considerable number of typographic errors in the text 

We use the notation and results of jl] 



The authors would like to thank Eamonn O'Brien and Guntcr Malle who found an error in the 
original version of the algorithm (Lemma 2.4 of M was false). I am also grateful to Sukru Yalcinkaya 



1.1 Description of the algorithm: the easy case, n is odd. 

Recall that ^l2n+i{q) denotes the kernel of the spinor norm 

S02„+i(g) ^{±1}; 

it is a subgroup of index 2 in S02n+i{q) and coincides with the commutator of S02n+i{q) ii n ^ 2 and 
(7^3 (for details see, for example, Artin j^). It is well known that PSp2((7) — ^siq) — PSL2((j) and 
PSp4(g) ~ ^siq)- Hence it is obviously enough to give an algorithm for n ^ 3. 

The easy case of the algorithm deals with the situation when n is odd, and the hard case is concerned 
with even n. 

If X is an element of even order, we use notation i(a;) for the involution in the cyclic group (x). 
The structure of the algorithm in the easy case is the same as described in [Q and very simple. 

* Research supported by the Royal Society and Lever hulme Trust. 
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INPUT: A black box group G isomorphic to PSp2„((j') or 172„+i(g), n ^ 3 and q > 3 given, n is odd. 

OUTPUT: If G is PSp2„((7) then, with probabihty > 1/2, the algorithm produces a witness against 
isomorphism G ~ f^2n+i(9)- 

(1) Find an element x of good order (the definition is given below). As it was discussed in ||l[], this can 
be done in polynomial time. 

(2) Compute the involution i = i(x) in the cyclic group generated by x, this also can be done in 
polynomial time. 

(3) This step is vacuous in the easy case. 

(4) Compute y :— i^i'^ for random elements g,h € G. If o{y) divides q{q + 1) or q{q — 1) then go 
to Step 5. Otherwise return 'G is PSp2„(g) ' , with the pair {i^,i'^) being a witness against the 
isomorphism G ~ Q2n+i{Q)- 

(5) Repeat Step 4 at least 20n times. If o(y) always divides q{q + 1) or q{q — 1) then return 'G is not 
PSp2„((7) '. The error probability is < 1/2 in this case. 

The number of necessary repetitions in Step 5 is a very crude estimate, which was sufficient for our 
theoretical purposes. In practice, if G is of type C„ and q > 3, then one does get an element i^ ■ i^ whose 
order does not divide q{q + 1) or q{q — 1) at the very first iteration. 

1.2 Description of the algorithm: the hard case, n is even. 

Let p be a prime number q = and m ^ 2. A prime dividing q™ — 1 but not — 1 for 1 ^ i < m is 
called a primitive prime divisor of q"^ — 1. It will be convenient to say that p is a primitive prime divisor 



of + 1 if p is a primitive prime divisor of g^™ — 1, see Lemma2.2 for further discussion of primitive 
prime divisors. In our case primitive prime divisors always exist by Zsigmondy Q. The primitive divisor 
rank pdrank(f;) is the maximal integer k such that o{g) is divisible by a primitive prime divisor of q^ — 1. 

INPUT: A black box group G isomorphic to PSp2„(g) or f22n+i(<z), n ^ 3 and q given, n is even (and 
hence q" = 1 mod 4) . 

OUTPUT: If G is PSp2„((j') then, with probability > 1/2, the algorithm constructs a witness against 
the isomorphism G ~ ^2n+i{Q)- 

(1) Find an element x of good order as in the easy case. 

(2) Compute the involution i — \(x) in the cyclic group generated by x. 



(3) Our next task is to make sure that i is an involution of type ii, i„_i or See Table y for the list 
of conjugacy classes of involutions and their centralisers (the details are mostly taken from [^). 

For this purpose, construct (almost) uniformly distributed random elements of Gdi) as explained 
in (l| Section 3]. The possibilities for Gcii) are given in Table |l|. Notice that the groups Caiti), 
Ccitn-i) (in the group G — ^l2n+i{q)) and Gcitn) have much bigger orders of semisimple elements 
than centralisers of involutions of other types. Namely, in these three cases we can find in Cg(i), 
with probability > 0{l/n), a 'big' element, that is, an element of order divisible by a prime 
primitive divisor of g"^^±l. However 'big' elements are not present in the centralisers of involutions 
of other types. As soon as a 'big' element in Cc{i) is found, we conclude that i is of type ti, tn-i 
(in the group G = r^2n+i(9)) or t„ and move to Step 4. If we cannot find a 'big' element, we return 
to Step 1 and start our search for an involution of type ti, or t„ from the beginning. 
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Involution ti \ Group G 


PSP2„(9) 


^^2n+l(9) 


h 


Sp2(g) oSp2„_2(?) 


Zq±i X f22n-l(g) 


ti 

2 sc: z sC n/2 (in PSp„((z)) 
2 < i < n - 1 (in ^l2n+i{q)) 


SP2»(9)°Sp2„_2,(g) 


n+^{q) X f^2(n-i)+i(g) 

(i even) 

f]~(g) X f^2(n-i) + l('7) 

(« odd and = —1 mod 4) 


tn-l (in Cl2n+l{q)) 




n2+„-2(9) X PSL3(g) 
(n is odd) 

(n is even and = — 1 mod 4) 


q = 1 mod 4 






g = — 1 mod 4 


^GU„((?2) 





Table 1: Involutions in groups PSp2„((j') and fl2n{<l), n even. The second and the third column show the 
generalised Fitting subgroup F* {Ca(ti)) of Coiti) in the groups PSp2„(g) and il,2n+i{(i), correspondingly. 
In all cases CciU) / F* {Cciti)) is an elementary abelian group of order < 4. 
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(4) Now we want to make sure that we found an involution of type t„. For that we compute pdrank(Mf ) 
for 20n random elements g. We know from Lemmata ^^8, 2.1C and 2.12 that if i is not of type t 



or if i is of type i„ and G — fl2n+ii<l) then pdrank(M^) ^ 6. Hence if we found g G G such that 
pdia.nk{ii^) > 6 then we have a definite witness against G being the orthogonal group. In this case 
we return ' G is Sp2„(q)'. 

If we have not found an element g with pdrank(M3) > 6, we can return to the search from Step 1. 

(4bis) Shortcuts. However, there are some shortcuts which give a quicker, although not one-sided 
solution. 



(4i) It follows from Lemma p.l2| that if pdrank(M3) ^ 6 for 0{n) random elements g & G and, for 
some g, pdrank(ii^) > 4 then i is likely to be an involution of type tn-i in the orthogonal 
group ri2n+i(<7)- 



(4ii) It follows from Lemmata 2.10 and 2.12 that if pdrank(M3) ^ 4 for 0{n) random elements 
g G G and, for some g, pdrank(M^) > 2 then i is likely to be an involution of type ti in the 
orthogonal group fl2n+i{q)- 



(4iii) It follows from follows from Lemmata ^.IC and 2.12 that if 



pdrank(M^) ^ 2 for 0{n) random elements g & G, then i is likely to be an involution of 
type ti in Sp2„(g') or of type t„ in fl2n+iiq)- Since the centralisers of these involutions have 
very different structure (see Table |l|), a further distinction between the two cases can be made 
by running a black box for C — Ccii)- The following steps allow to detect a normal subgroup 
SL2((7) in C = Csp2^(q)(ti). 

* Choose e — ±1 so that — e = 2 mod 4. Randomly search in C for an element x 
of order divisible by 4 and by a primitive prime divisor of — e. If such an element 
cannot be found after 0{n) attempts, i is likely to be an involution of type i„ in ^l2n+i{<l)- 

* Once the element x is found, compute an element y of order 4 in {x). 

* Compute 0{n) elements yy'^ for random c e C. If, for all of them, o{yy'^) \ q{q±l) then i is 
likely to be an involution of type ti in T'Sp2nil) (a-nd y belongs to the SL2((j')-coniponent 
of Cg«). 

* If, however, o(yy'^) \ q{q ± 1) for some c £ C then i is likely to be an involution of type t„ 
in ^2n+iiq)- 

(5) If we returned to Step 1 more than 20n times and have not found a witness against the isomorphism 
G ~ Q,2n+iiq), we conclude that G ~ rt2n+i{q)- 

2 Tori and semisimple elements in groups and C„ 

Definitions and notation used in the paper are standard and borrowed from Carter ||^, Q and Gorenstein, 
Lyons and Solomon [||. We retain the notation and use the results of [Q. We use the description of 
maximal tori in groups i?„ and C„ as given in (l[ Section 2.1]. 

2.1 Tori of order ± 1). 

The crucial steps of our algorithm are aimed at finding involutions belonging to tori ^Tqn±i twisted by 
an element of the cycle type ±n. We have to distinguish the following cases: 

(i) q" = —1 mod 4, i.e. n odd and q = —1 mod 4. 

(m) q" = 1 mod 4, i.e. n even or q = 1 mod 4. 
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In the first case the maximal tori of order ■^{q"' + 1) contain involutions, in the second case the maximal 
tori of order ^(g" — 1) contain involutions. 
We call tori ■^Tqn±i maximal twisted tori. 

We are looking for involutions belonging to maximal twisted tori. For that purpose, we introduce 
the concept of 2-height. Let X be a group and x ^ X such that o{x) = 2*''m where m is odd. Then k is 
called the 2-height of x. The 2-height of X is the maximum of heights of elements in X. 

Lemma 2.1 If T is a torus of type iTgii±i then 

iNciT) : T\ = 2n. 



Proof See [0, Lemma 2.3]. □ 

2.2 Primitive prime divisors. 

Let p be a prime number p and m ^ 2. A prime dividing p™ — 1 but not — 1 for 1 ^ i < m is called 
a primitive prime divisor of p'^ — 1. In our case primitive prime divisors always exist by Zsigmondy |^]. 
Furthermore an integer j is called a ppd{p] m) if j \ q™ — 1 and j is divisible by a primitive prime divisor 
of — 1. Let q = p'^. A primitive prime divisor or p"^ — 1 will be called, for brevity, a ppd{q, m). An 
element of G of order dividing a ppd{q, n) is called a primitive prime divisor element. 

The role of primitive prime divisors becomes apparent from the following simple lemma. 

Lemma 2.2 (i) Let r be a ppd{q,n). If r \ q™ + 1 for m < n, then n is even and m = n/2. 

(ii) Let r be a ppd{q, 2n). Then r \ ± 1 for all m < n and r j — 1. 

Proof, (i) Assume that r | + 1 for some m < n. This means that r | g^" — 1 and the Galois field 
F^2m has a multiplicative subgroup R of order R. But the Galois field ¥gn also contains a multiplicative 
subgroup of order r. Since the algebraic closure contains a unique field of every give order, the 
intersection Fq2m n Fgn is well defined. The inclusion ¥q,^ < ¥g2m because this would mean n \ 2m, 
contradicting our assumption m < n. The equality ¥g2m — ¥qn leads to the desired conclusion m = n/2. 
Finally, we are left with the case when F^fc = ¥g2m n Fg^ is a proper subfield of Fgn . But then the 
multiplicative group of F^k contains _R, an unique subgroup of order r in the both Fg2m and ¥qn , and 
therefore r | q'^ — 1, a contradiction. 

(ii) follows from (i); we need only notice that r is odd, and, since q" — 1 and g" + 1 have no odd 
divisors in common, r | q" — 1 □ 

2.3 Good elements. 

Let P be a maximal twisted torus of even order and let \P\ = 2*^771 where k and m are natural numbers 
and m is odd. We say that an element 5 € P is of maximal 2-height, if the 2-height of g is fc. Since 
P is cyclic, elements of maximal 2-height exists. Furthermore half of the elements of P have maximal 
2-height as the cyclic group generated by a 2-element x of maximal 2-height contains 2'^"^ elements of 
maximal 2-height, namely a;" for s an odd integer. 

In the following we are looking for elements of maximal 2-height in a maximal twisted torus. We call 
g good if it satisfies the following conditions: 

• pdrank(g) = 2n when q" = — 1 mod 4 and pdrank((7) = n when q" = 1 mod 4. 

• g has maximal 2-height. 

We shall call good elements belonging to maximal twisted tori in G good twisted elements. 
Lemma 2.3 The number of good twisted elements in G is at least -g^lGI. 
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Proof. Let T be a maximal twisted torus in G. By Lemma 2.1 |A''G'(r)/T| = 2n. Therefore there are 
2^q^ tori conjugate to T. Let r be a ppd{q, n) or ppd{q, 2n), depending on whether T is of type ^Tq^^i 

or iTgn+i, and denote by T^^ the unique subgroup of index 2r in T. Then any element whose image 
in T/T"^^ generate T/T"^^ is a good element. It is easy to see that the number of such elements in T is 
^^|r|. Also, it follows from Lemma 2.2 that a good element does not belong to a proper subtorus of T 
and thus belongs to a unique twisted torus. Now we can estimate the total number of good elements in 
all twisted tori as at least 



2nT 



2r 



An 



5n 



since, obviously, r ^ 5. 



□ 



2.4 Regular twisted elements. 

A regular element in a torus of type iT^,i±i or will be called a regular twisted element. Notice that good 
elements in these tori are regular twisted elements. 

Lemma 2.4 The number of regular twisted elements in G conjugate to an element in a torus T of the 
given type ^Tqn_i or ^Tgn+i is at least |^|G|- 

Proof: See |], Lemma 2.5]. □ 

2.5 Involutions of type t„ in PSp2„(<?). 

For the following compare the description of conjugacy classes of involutions in PSp2„(fc) and S02,i+i(fc) 
as given m Theorems 11.52 and 11.53 of Taylor H and Section 4.5 of |||. 

The group PSp2„ (q) has [^] + 1 conjugacy classes of involutions. Let s G Sp2„ (q) represent an element 
of order 2 in PSp2„((i'). Then either — ~I or ~ I. All involutions which are represented by an 
element s such that = — / form one conjugacy class of involutions of type /;„. Furthermore for every 
even number 2 ^ m ^ n there is one conjugacy class of involutions of type tjn/2- If s is an involution of 
type tm/2j then the eigenvalues of s on the underlying natural module are —1 with multiplicity m and 
1 with multiplicity 2n — m. 

Lemma 2.5 Let G := PSp2„((j') and i Cz G an involution of type If H is a torus of G then there 
exists an involution j £ i^ such that H is inverted by j . 



Proof See || Lemma 2.6]. □ 

Lemma 2.6 Let G :~ PSp2„(9) and i G G be an involution of type tn- Then i^i^ is a regular twisted 
element of odd order with probability p ^ 2^ /c" random g,h G G. 

Proof See [§ Lemma 2.7]. □ 
As an immediate corollary we have the following lemma. 

Lemma 2.7 If i is an involution of type tn in G — PSp2„((z) then the share of elements g £ G such 
that the element ii^ is of odd order is at least Trk- ■ 

20n 

Proof: See |l Lemma 2.8]. □ 
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2.6 Involutions of type t„ in ^^2n+i{(l)- 

The group ^2n+i{Q) has n conjugacy classes of involutions. There is one conjugacy class of involutions 
of type t„j/2 for each even number 2 ^ m ^ 2n. An involution s is of type t„i if the eigenvalues of s on 
the underlying natural module are —1 with multiplicity m and 1 with multiplicity 2n + 1 — m. 

Lemma 2.8 Let G :— il2n+i('z) o-nd i ^ G an involution of type i„. Then the element i^ ■ i^ acts 
trivially on a subspace of codimension 2 in the natural module for G. In particular, o{i^i^) \ q{q + 1) or 
o{iH'^) I q{q - 1) for all g, heG. 

Proof. Let V be the natural module for G. Since i is of type t„ the dimension of the eigenspace V- of 
i for the eigenvalue —1 is 2n, while the eigenspace V+ for the eigenvalue +1 is 1-dimensional. Therefore 
i^i'' belongs to C = CsL(y)(VL(7 fl V-h) and W — V-g DV-h has codimension ^ 2 in V^. Notice that 
C/Op{C) is a subgroup of SL2(g) (here p = charF^). Notice also that the element i^i'^ leaves invariant 
the subspace V+g + V+h. The result now easily follows. □ 

Lemma 2.9 Let G :~ f^2n+i(<z) o,rid i £ G an involution of type tn- Then the share of elements g (z G 
such that the element ii^ is of odd order is at least 

1 

4(l + 2/g + o(l/g))- 



Proof: We use notation of the previous lemma with W — HVLg and U — V+ + V+g. Let q = 1 mod 4. 
Then a direct computation with the help of information about finite orthogonal geometries contained in 
Section in.6 of Artin shows that, with probability 

1 

P> 



4(l + 2/g + o(l/g))' 

the subspace U is nondegenerate and has index 0. Hence the group {i,i^) is isomorphically imbedded 
into 0{U) ~ Zg+i XI Z2 and half of the elements ii^ have odd order, yielding the desired estimate. The 
case q = —1 mod 4 can be treated analogously. □ 

2.7 Involutions of type ti in PSp2„(q') 
Lemma 2.10 // i is an involution of type ti in PSp2„((j') then 

o{iH'') I q{q±l). 



Proof: Similarly to the previous section, a preimage in G = Sp2„((7) of the involution i (which we denote 
by the same symbol i) is a 'low dimensional' involution, that is, the eigenspace V- has dimension 2 or 
n — 2. Assume that dimy_ = 2, the other case is similar. Then, for random g,h G G, the group {i^,i^) 
leaves the subspaces V^g + V-h and V+g fl V+h invariant and V = {V-g + V-h) © {V+g n V+h). Denote 
W = V-g + V-h. Notice that {is,i'') < Sp(Vr). 

We have three possibilities: dimVF = 2, 3 or 4. If diniM^ = 2 the i^ — If the restric- 
tion of the symplectic form to W is degenerate (this automatically holds when dimVF = 3) then 
Sp(M^)/Op(Sp(W^)) ~ Sp2((7) = SL2(g), and the lemma follows. If dimVF = 4 and the restriction 
of the symplectic form to W is non-degenerate then Sp(H^) = Sp4(q) and i^ and i^ are involutions of 
type ti in Sp(VK). Now the result follows from an easy computation in Sp4(<;). □ 

Moreover, the same analysis yields the following result. 

Lemma 2.11 // i is an involution of type ti in PSp2nil) then the .share of elements in G .such that ii^ 
is of odd order is at least 1/4. 
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2.8 Involutions of type ti and in ^l2n+i{(l) 

Lemma 2.12 Let i be an involution of type ti or in G = ^2n+i{Q) o-nd g is a random element in 
G. 

(i) // i is of type ti then pdrank(ii5) ^ 4. Moreover, pdrank(M5) > 2 with probability bounded from 
below by a constant. 

(ii) // i is of type i„_i then pdrank(M^) ^ 6. Moreover, pdrank(M3) > 2 with probability bounded from 
below by a constant. 

(iii) // i is of type ti or tn-i then the element ii^ has odd order with probability bounded from below by 
a constant. 

The experimental evidence suggests that the constants in (i) and (ii) are 1/8 for ti and 1/4 for tn-i. 
Proof. The proof is similar to the previous result. □ 

2.9 Involutions produced by good elements. 

Lemma 2.13 Let i be the involution in a maximal twisted torus T < G . Then i is of type tn. 

Proof See H Lemma 2.13(1)]. □ 



We shall use the following corollary from this result and Lemma 2.3 



Lemma 2.14 Let g be a random element in G. Then the involution i produced by g has type tn with 
probability at least l/5n. 

Finally, 

3 Detailed description of the algorithm 

3.1 Step 1: search for a good element x. 

The algorithm spends most of its timing testing if elements are good. This is done exactly as described 
m 0, Section 4]. 



It easily follows from Lemma 2.3 that a sample of N random elements of G contains a good element 
with probability 



1 - ( 1 - — I >l-e-- 

Therefore if we are given fixed margin of error e, < e < 1, and take a sample of size 0(7ilog(l/e)), then 
we can be satisfied that it contains a good element with probability 1 — e. As proven in [|l], production 
of a random element and testing it for goodness require polynomial time. 

3.2 Step 2: involution i = i{x). 

As soon as we have a good element x and know that it is an element of maximal 2-height in a cyclic 
group of order ^{q'^ i l)i the involution x = i(a;) can be o btain ed by raising x to the power ^(q" ± 1). 
0, 



As discussed in 111, this requires polynomial time. Lemma 2.13 asserts that i is of type i„ 
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3.3 Step 3: involution of type 

Notice that if n is odd then the involution i is of type <„. 

However, if n is even then i can happen to be of different type. 

We can check whether i is of type i„ or not by constructing uniformly distributed random elements 
Ciix) of Ccii) as described in Section 3]. If z is of type ii, t„_i (in fl^j^iq)) or f„ then, in view of 
Lemmata ^.11 , 2.12 and 2.£, this method produces (in polynomial of nlogg time) random and uniformly 
distributed elements of Ca{i). Let L — F * {CG{i)). Depending on whether i is one of the ti, t^-i (in 
^tnil)) or tn, L contains 'big' elements, that is, elements of order divisible by a primitive prime divisor 
of - 1, + 1, + 1 (compare Table |l] and Table |). It can be observed directly from the 
formulae for orders in Table ^jthat centralisers of involutions of other types do not contain 'big' elements. 



\GU{q)\ 




|PSp2fe('?)l = \^2k+i{q)\ 


(5^g'=^(g2_l)(g4_i)...(g2fe_i) 




_l_^Mfc-l)(q2 _ i)(^4 _!)... (g2fe-2 _ _ 1) 




_l_qk(k-l)^q2 _ i)(^4 _!)... (g2fc-2 _ + 1) 


|GUfc(q2)| 


^fe(fe-l)/2(^ + - l)(q^ + l)(g4 _ 1) . . . (^fe _ 



Table 2: Orders of some classical groups. 



pp. 



It is easy to see that 'big' elements are regular elements in maximal tori of orders ^((? ± 1)(<7"^^ ± 1) 
or \[q"' + 1). Using the same technique as in Sections 2.1-2.3, we can see that the probability of a 
random element to be 'big' is 0(l/n) and hence one of them has to pop up from the black box for Ccii) 
after 0{n) iterations. As soon as we have found a 'big' element, we know that i is of type ti, (in 
^2n+i{q)) or tn and can go to Step 4. 

If a 'big' element cannot be found after 0{n) tests of elements from Ccii), then to Step 1 and start 
the search for an involution of type ti or i„ from the beginning. 



3.4 Step 4. 

Now we want to make sure that we found an involution of type t„ in PSp2„(q). For that we compute 
pdrank(M^) for 20n random elements g. We know from Lemma 2.6 that if i is of type tn in PSp2„(q), 
then, with probabilit y > 1 — e ~^, t he pr imitive divisor rank of one of the elements ii^ is at least 2n. We 
know from Lemmata 2.8, 2.1C and ^.12 that if i is not of type tn, or if i is of type tn and G = ^}2n+i{q) 
then pdrank(M3) ^ 6. Hence if we found g & G such that pdrank(M^) > 6 then we have a definite 
witness against G being the orthogonal group. In this case we return ' G is Sp2nilY- 

If we have not found an element g with pdrank(M3) > 6, we can return to the search from Step 1. 

The estimate of the number of elements ii^ needed in this step is, of course, very crude; an experi- 
mental evidence suggests that, instead of 20n elements ii^, 8 elements suffice. 

Shortcuts in Step 4 are mostly self-explanatory. 
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3.5 Step 5. 

The involution i of type i„ behaves very differently in PSp2„(q) and Q2n+iiQ)- Namely, in PSp ^„(g ) 
an element i^i'^ for random g,h € G is a. regular twisted element with probability ^ (Lemma 2.6), 
while in U2n+i{q) we always have o{i^i^) \ q{q ± 1) (Lemma |2.8|) . Obviously, a regular twisted element 
cannot have order dividing q{q ± 1). As soon as we have that o(i^i^) does not divide q{q ± 1), the pair 
(i^ , i^) is a witness against the isomorphism G ~ n2n+i{q) and in that case G ~ PSp2n{l)- 

On the other hand, if in 20n tests we had the same result o[i^i^) \ q[q ± 1), then G ~ ^2n+i[q) 
with probability ^ 1 — e^^. Timing estimate here is very crude and can be considerably improved. In 
practice, if G ~ PSp2„((7) then o{i^i^) \ q{q ± 1) at the very first attempt. 



4 The GAP code. 

The algorithm was tested in GAPQ. The Appendix contains a (very primitive) GAP source code Bn-v-Cn . g, 
which generates (at random) a group G = Sp2nil) '^^ S02n+i('?) and initiates the functions Test(G) 
and PlainTest (G) which determine the type of the group G with the help of the algorithm described 
in this paper, with (Test(G)) or without (PlainTest (G) ) shortcuts in Step 4. 

The functions PlainTest (G) and Test (G) return the list of 4 integers. Of these, 

• The first position is or 1 , according to the answer symplectic or * orthogonal given by the algorithm. 



The second position is the rule number used in this decision (used for the debugging of the algo- 
rithm) . 



• The third position is the dimension of the eigenspace for the eigenvalue —1 of the involution used 
in the identification of the group this information is used only for the debugging of the algorithm) . 

• The fourth number is the numbers of rc-runs of the test if the first run was inconclusive. 

I wish to emphasise that this GAP code was written only for checking the algorithm and is not 
designed for the practical use. In particular, it works with the groups Sp2„(g) and S02n+i('7) instead of 
PSp2„((7) and f22n+i('7), although all computations are done in the way which makes sense in the factor 
group modulo the scalars and in the commutator group; in particular, it works, instead of involutions, 
with pseudoinvolutions, that is, elements of order 2 or 4 whose images modulo the scalars are involutions. 

To see how the code works, enter at the GAP command prompt: 

gap> Read("Bn-v-Cn.g") ; test : =Test (G) ; Interpret (test) ; 



gap> Read("Bn-v-Cn.g") ; test : =PlainTest (G) ; Interpret (test) ; 

To check the answer, enter 
gap> G; 

I have to warn, however, that PlainTest (G) takes considerable time to produce the answer if G is 
the orthogonal group. 

So far more than 2000 experiments with the function Test (G) have not yield a single wrong answer, 
while PlainTest (G) erred once in 600 experiments. 

This was my first ever attempt to write in GAP, and I bravely ignored the good programming practice; 
in particular, the global variables q, n, type (0 or 1 for symplcctic or orthogonal) for generation of a 
black box group can be changed only in the file Bn-v-Cn.g, which after that has to be read again to 
generate all the required functions. 



^The software package GAP is available, for variety of platforms, from ittp://aldebaran.math.rwth-aachcn.de~GAP/ 
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Also, I am using the function Order (element) , to make the code simpler; this is not necessary, of 
course, and it can be replaced by raising the element in the appropriate powers as discussed in [|l|. This 
would make it possible to work in "large" groups. 

Structurally, I was not concerned with the optimal branching in the program. I intentionally wanted 
all branches being regularly visited, although, in many cases, this could be avoided. 
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Appendix: GAP code 



#################################################################### 
### 

### Bn-v-Cn.g 
### 

### A. V. Borovik, April 2001 
### 

### The algorithm tests involutions in orthogonal and symplectic groups 
### thus distinguishing between the two. 

### Given a linear group G = PSp(2n,q) or Omega(2n+l, q) , 
### with known q and n, 

### the algorithm determines whether G is sjrmpelctic or 
### orthogonal . 

### The algorithm implemented here contains a fix for a bug 
### in the original theoretic version, 

### C. Altseimer and A. Borovik, "Probabilistic recognition 

### of orthogonal and symplectic groups", 

### Groups and Computations III, (eds W.'M.Kantor and A. 

### Seress), de Gruyter, Berlin and New York, 2001, pp. 1-20. 

### 

### This code was not designed for the practical use, 

### but just for checking the theoretical algorithm. 

### In pEirticular, since I had no means for generation of 

### black boxes for PSp(2n,q) for large n and q, 

### the algorithm actually works with Sp(2n,q) instead of 

### PSp(2n,q), but make only those group calculations 

### which make sense in its factor group PSp(2n,q) . 

### In pairticular, it works with "pseudoinvolutions" , 

### elements of order 2 or 4 which have order 2 in G/Z(G) . 

### 

### Recommended use: 

### gap> Read("Bn-v-Cn.g") ; test :=Test (G) ; Interpret (test) ; 
### or 

### gap> Read("Bn-v-Cn.g") ; test :=PlainTest(G) ; Interpret (test) ; 
### To check the answer, enter 
### gap> G; 

#################################################################### 



### Global constants: q, n, type, G, effort 
q:= Random([5,7,9,ll,13,17,19,23,25,27,31]) ; 



if q >19 then 

n : = Random ( [3 . . 6] ) ; 

fi; 

if q > 10 and q < 20 then 
n:= Random( [3. . 12] ) ; 



fi; 

if 

fi; 



q < 10 then 
n:= Random( [3. .20] ) ; 



### selection of the dimension of 
### of the group with the aim to 
### keep the exponent at 
### computationally feasible size 
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type:= Random ( [0, 1] ) ; ### type = means symplectic 

### type = 1 means orthogonal 

if type = 1 then 
G:=S0(2*n+l,q) ; 

else 

G:=Sp(2*n,q) ; 

fi; 



### Generation of the group G, type =1 
### Generation of the group G, type = 



effort := 2; 



### This parameter determines how hard 
### we are looking for the witness 
### elements. Most seeirches make 
### effort * n attempts to find a 
### "big" element 

### Effort =4 corresponds to the theoretical threshold 
### in the group Sp(2n,q) or Omega ( 2n+ 1 ,q) , 
### a random element is a regular element 

### in a maximal twisted torus of order 

### q~n -1 or q"n+l with probability > l/4n. 



### 

### Functions 
### 



OddPart := function(nn) 

### Odd part of the integer nn 

local m; 

m :=nn; 

while m mod 2 = do 
m := m/2; 

od; 

return m; 

end; 



EvenPart := function(nn) ; 

### Even part of the integer nn 

return nn/OddPeirt (nn) ; 

end; 



PrimitiveDivisor := fmiction(qq,nn) 
### Finds the largest divisor of the number qq'nn -1 
### which has no non-trivial divisors in common with 
### any of the numbers qq~l - 1 for 1 < nn. 

local pd,ll ,gcd, ii; 

pd := qq~nn-l; 

for ii in [l..nn-l] do 

11 := qq~ii -1; 

gcd := Gcd(pd,ll) ; 

while gcd > 1 do 
pd := pd/gcd; 
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gcd:= Gcd(pd,gcd) ; 

od; 

od; 

return pd; 

end; 



Pseudolnvolution := function(g) 

### Finds an element in the cyclic group <g> 

### which is am involution in G/Z(G) 

### The function uses the standaird GAP function Order (g) 
### which uses prime factorisation and does not work in big 
### groups. However, it is easy to eliminate calls to Order 
### and replace them by exponentiation. 

local order , i , rank ; 
order := Order (g) ; 
if order mod 2=1 then 
return g~0; 

else 

i := g* (order/2) ; 

fi; 

rank := REmk(i+i~2) ; 

if (rank = and type = 0) then 

if order mod 4 <> then 
return g~0; 

else i := g" (order/4); 

fi; 

fi; 

return i; 

end; 



InCentraliser := function(j) 

### Returns an element from the centraliser C_G(j) 

### Most of the elements produced by this function are 

### involutions; their distribution in C_G(j) is 

### invariant vmder action of C_G(j) by conjugation 

### Some (relatively small) portion of elements is 

### produced by a different algorithm $\zeta_l$ in 

### notation of our paper; these elements are distributed 

### imiformly in the group C_G(j). 

local gg,y,o,z,c; 

gg := PseudoRandom(G) ; 

y := j * j"gg; 

o := Order (y) ; 

if j~4 <> j~2 then 

if (o = 1 or o = 2 or o = 4) then 

z := y; 

fi; 

if o mod 4=0 then 

if y"(o/2) = j'2 then 
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z := y-(o/4); 

else 

z:= y-(o/2); 

fi; 

fi; 

if o mod 4=2 then 

if y~(o/2) = j~2 then 

z := gg*y-(((o/2)-l)/2); 

else 

z:=j; 

fi; 

fi; 

if o mod 2=1 then 

z := gg*y"((o-l)/2); 

fi; 

else 

if o mod 2=0 then 
z := y~(o/2) ; 

else 

z := gg*y'((o-l)/2) ; 

fi; 

fi; 

return z; 

end; 



OccurrencesOf PrimitiveDivisors := function (m) 

### Finds all integers k such that m is divisible by 

### a primitive prime divisor of q"k-l . 

local tuple, k; 
tuple := 0*[1. .2*n] ; 
for k in [1. .2*n] do 

if Gcd(PrimitiveDivisor (q,k) ,m) > 1 then 
tuple [k] : =k ; 

fi; 

od; 

return tuple; 

end; 

PrimitiveDivisorRankOf Element := function(g) 
### Finds the maximal integer k =< 2*n such that 
### the order of the element g is divisible 
### by a primitive prime divisor of q~n-l 
### (remember, q and n are global constants) . 
### Again, calls to Order (g) can be eliminated 

return Maximum (OccurrencesOf PrimitiveDi visors (Order (g) ) ) 

end; 



IsGood := function(g) 

### Tests whether the element g is good or not. 
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if (q"n -1) mod 4=0 then 

return (PrimitiveDivisorRankOf Element (g) = n and EvenPart (Order (g) ) = 

(EvenPart(q"n-l)/(2"type))) ; 

fi; 

if (q"n + 1) mod 4=0 then 

return (PrimitiveDivisorRankOf Element (g) = 2*n and EvenPeirt (Order (g) ) = 

(EvenPart (q"n+l) / (2"type) ) ) ; 

fi; 

end; 

GoodElement := function(G) 

### Finds a good element in the group G, 

### of course, under the assumption that 

### G = Sp(2*n,q) or S0(0,2*n+l,q) ; 

### in the latter case, the good element 

### is actually found in the commutator group 

### 0mega(2*n+l,q) of S0(0,2*n+l,q) . 

local gl, g2, g; 
if type = 1 then 

gl : =PseudoRandom(G) ; 
g2:=PseudoRandom(G) ; 
g :=gl--l*g2--l*gl*g2; 
while IsGood(g) = false do 
gl : =PseudoRandom(G) ; 
g2:=PseudoRandom(G) ; 
g :=gl--l*g2--l*gl*g2; 

od; 

else 

g:=PseudoRandom(G) ; 
while IsGood(g) = false do 
g:=PseudoRandom(G) ; 

od; 

fi; 

return g; 

end; 



TypeOf Centraliser := function(j) 

### Returns one of the lists [] , [1], [2], and [1,2], 

### depending on whether an element of order dividing 

### a primitive prime divisor of q~(n-l)-l or a 

### primitive prime divisor of q~(n-l)+l is found in C_G(j) 

### (this adds 1 or 2 to the list , correspondingly) . 

### 

### The procedure is pretty slow; it can be made faster 

### by a better way of generating random elements in C_G(j) 
### (for example, the use of the Product Replacement Algorithm 
### should considerably speed it up) . 

local X, maxppd, typeOfC, counter; 
X := InCentraliser(j) ; 
maxppd : =0 ; 
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typeOfC := Set([]) ; 
counter := 0; 

while typeOfC <> Set ([1,2]) and counter < effort*n do 

X := X * InCentraliser ( j ) ; 

maxppd := PrimitiveDivisorRankOf Element (x) ; 
if maxppd = n-1 then 
AddCtypeOfC, 1); 

fi; 

if maxppd = 2* (n-1) then 
Add(typeOfC, 2); 

fi; 

counter := coimter +1; 

od; 

return Set (typeOf C) ; 

end; 

ProductsOf Conjugates := function(j) 
### Determines how big is the primitive divisor rank 
### of the product j * j~x of the element j and 
### its random conjugate j~x. 

local X, x2, y, maxppd, coimter; 
maxppd := 0; 
counter := 0; 
if type = then 

while (maxppd < 8 and counter < effort*n) do 

x := PseudoRandom(G) ; 

y := j * j-x; 

maxppd := Maximum ( [maxppd, PrimitiveDivisorRankOfElement(y)] ) ; 
coimter := counter +1; 

od; 

else 

while (maxppd < 8 and counter < effort*n) do 
X := PseudoRandom(G) ; 
x2 := PseudoRandom(G) ; 
X := x*-l * x2~-l ♦ X * x2; 

y := j * j-x; 

maxppd := Maximum ( [maxppd , PrimitiveDivisorRankOfElement(y)] ) ; 
coimter := counter +1; 

od; 

fi; 

return maxppd; 

end; 



Criticallnvolution := function(G) 

### Finds an involution j in the cyclic group generated by a good element, 
### with the extra property that the centraliser C_G(j) has an element 
### whose order is divisible by a primitive prime divisor of q*(n-l)-l 
### or q-(n-l)+l. 

local X, j, typeOfC; 

typeOfC := []; 

X : = GoodElement (G) ; 

j := Pseudolnvolution(x) ; 

while (typeOfC = Set([]) or j=j'2) do 
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X := GoodElement (G) ; 

j := Pseudolnvolution(x) ; 

typeOfC := TypeOfCentraliser(j) ; 

od; 

return j ; 

end; 



PreTest := function(G) 

### Tests the group G for being symplectic or orthogonal 
### The result may be inconclusive, in which case 
### the test has to be repeated. 

local g, j, typeOfC, conj , conj2, x, y, z, minus, plus, majxppd, 
evenpeirt, counter, maxppd2 , coimter2 ; 

### The case when n =3 

if n =3 then 

g : = GoodElement (G) ; 
j := Pseudolnvolution(g) ; 
conj:= ProductsOf Conjugates (j ) ; 
counter :=0; 

while conj < 3 Eind counter < effort * n do 
conj2 := ProductsOf Conjugates(j) ; 

conj := Maximum ( [conj , conj2]); 
counter := coiuiter + 1; 

od; 

if conj < 3 then 

return [1, 13, Rank(j-j"2)/2] ; 

else 

return [0, 13, Rank(j-j"2)/2] ; 

fi; 

fi; 

### The general case 
j := Criticallnvolution(G) ; 
conj := ProductsOf Conjugates(j) ; 
if conj > 7 then 

return [0, 1, Rank(j-j-2)/2] ; 

fi; 

typeOfC := TypeOfCentraliser(j) ; 
if typeOfC = [1,2] and conj = 4 then 
return [1,2, Rank(j-j"2)/2] ; 

fi; 

if n mod 4=2 and q mod 4 =1 and typeOfC = [1] and conj = 6 then 
return [1,3, Rank(j-j~2)/2] ; 

fi; 

if n mod 4=2 and q mod 4 =3 and typeOfC = [2] and conj = 6 then 
return [1,4, Rank(j-j*2)/2] ; 

fi; 

if n mod 4=0 and typeOfC = [1,2] and conj < 3 then 
return [1,5, Rank(j-j~2)/2] ; 

fi; 
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if n mod 2=1 and typeOfC = [1,2] and conj < 3 then 
return [1,6, Rank(j-j"2)/2] ; 

fi; 

if n mod 2=0 and typeOfC = [1,2] and conj < 3 then 
#### 

minus := (EvenPart(q~ (n-l)-l) mod 4 <> 0) ; 
plus := (EvenPart(q"(n-l)+i) mod 4 <> 0) ; 
X := InCentraliser(j) ; 

if minus then 
maxppd : =0 ; 
evenpart : =1 ; 

counter := 0; 

while (maxppd <> n-1 or evenpart < 4) and counter < 3*effort*n do 
X := X * InCentraliser(j) ; 
meixppd := PrimitiveDivisorRankOf Element (x) ; 
evenpeirt := EvenPart (Order (x) ) ; 
counter := counter +1; 

od; 

if (maxppd <> n-1 or evenpart < 4) then 
return [1,7, Rank(j-j*2)/2] ; 

fi; 

fi; 

if plus then 
maxppd : =0 ; 
evenpart : =1 ; 
counter := 0; 

while (maxppd <> 2* (n-1) or evenpart < 4) and counter < 3*effort*n 

X := X * InCentraliser ( j ) ; 

maxppd := PrimitiveDivisorRankOf Element (x) ; 
evenpart := EvenPart (Order (x)) ; 
counter := coimter +1; 

od; 

if (maxppd <> 2* (n-1) or evenpart < 4) then 
return [1,8, Rank(j-j"2)/2] ; 

fi; 

fi; 

y := x~ (Order (x)/evenpart) ; 

counter2:=0; 

maxppd2 : =0 ; 

while (maxppd2 < 3 and counter2 < effort*n) do 
X := X * InCentraliser (j) ; 

z := y * y"x; 

maxppd2 := Maximum ( [maxppd2 , PrimitiveDivisorRankOfElement(z)] ) ; 
co\mter2 := coiuiter2+l; 

od; 

if maxppd2 < 3 then 

return [0,9, Rank(j-j*2)/2] ; 

fi; 

if maxppd2 > 2 then 

return [1,10, Rank(j-j"2)/2] ; 



A. V. Borovik • Orthogonal and Symplectic Groups, Revisited • 21.10.01 



fi; 
#### 

fi; 

if n > 3 then 

return [-1,11, Rank(j-j"2)/2] ; 

fi; 

end; 

Test := fvmction(G) 

### PreTest(G) is repeated until we get a conclusive answer 

local pretest, counter; 
pretest := PreTest(G) ; 

if pretest [1] = or pretest [1] = 1 then 

Add (pretest , 0) ; 
return pretest; 

fi; 

counter : =0 ; 

while pretest [1] = -1 do 
pretest := PreTest(G) ; 
counter := counter + 1; 

od; 

#Add (pretest, 0); Add (pretest, 0); Add (pretest, 0); 
Add (pretest, counter); 
return pretest; 

end; 



PlainPreTest := function(G) 



### Tests the group G for being 

### symplectic or orthogonal. 

### No shortcuts in the algorithm. 

### The result may be inconclusive, 

### in which case the test has to be repeated. 



local g, j, typeOfC, conj ; 
g : = GoodElement (G) ; 
#g := PseudoRandom(G) ; 
j := Pseudolnvolution(g) ; 
while j = j ~0 do 

g : = GoodElement (G) ; 

#g := PseudoRandom(G) ; 

j := Pseudolnvolution(g) ; 

od; 

conj := ProductsOf Conjugates (j ) ; 
typeOfC := TypeQfCentraliser(j) ; 

if (conj > 7 and typeOfC <> []) or (conj > 2 and n = 3) then 
return [0, 14, Rank(j-j"2)/2] ; 

else 
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return [-1, 14, Rank(j-j*2)/2] ; 

fi; 

end; 



PlainTest := functioii(G) 



### PlainPreTest(G) is repeated until we get a conclusive answer. 



local pretest, counter; 
pretest := PlainPreTest(G) ; 



if pretest [1] = or pretest [1] = 1 then 
Add(pretest,0) ; 
return pretest; 

fi; 

counter : =0 ; 

while pretest [1] = -1 and counter < ef f ort*n do 

pretest := PlainPreTest (G) ; 
counter := counter + 1; 

od; 

#Add(pretest, 0); AddCpretest, 0); AddCpretest, 0) ; 
Add (pretest, counter); 
if counter > effort*n-l then 
pretest [1] := 1; 

fi; 

return pretest; 

end; 



Interpret := function(test) 

### Interpretation of the test result 



local 1; 

1 := Length(test) ; 
if test[l] = then 

PrintC The group is symplectic. \n" ) ; 

fi; 

if test[l] =1 then 

PrintC The group is orthogonal. \n" ) ; 

fi; 



if test [2] = 1 then 

PrintC This is the DEFINITE answer because the test has found an involution j \n"); 

PrintC whose centraliser C_G(j) contains an element of\n") ; 

PrintC order divisible by a primitive prime divisor of q~(n-l)-l \n"); 

PrintC or q~(n-l)+l, and such that, in addition, there is a product\n"); 

PrintC j*j*x of two involutions conjugate to j whose order is divisible\n") ; 

PrintC by a primitive prime divisor of q"k-i for some k>6;\n"); 

PrintC this is possible only for an involution of type t_n in the sjnnplectic group. \n"); 

f i; 



if test [2] = 2 then 

PrintC An involution j was found whose centraliser C_G(j) contains\n") ; 

PrintC elements of order divisible by primitive prime divisors of the both\n"); 
PrintC q'(n-l)-l or q"(n-l)+l, and with the maximum of the primitive prime\n"); 
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Print (" 
Print (" 
Print (" 
Print (" 



if test [2] 
Print (" 
Print (" 
Print (" 
Print (" 



if test [2] 
Print (" 
Print (" 
Print (" 
Print (" 



if test [2] 
Print (" 
Print (" 
Print (" 
Print (" 
Print (" 
Print (" 



if test [2] 
Print (" 
Print (" 
Print (" 
Print (" 
Print (" 
Print (" 
Print (" 



if test[l] 
Print (" 
Print (" 
Print (" 
Print (" 



if test[l] 
Print (" 
Print (" 
Print (" 
Print ( " 
Print (" 



divisor ranks of random products of conjugates j~x * j~y being \n"); 
equal 4. This almost certainly identifies this involution j \n"); 
as involution of type t_l in S0(2n+1 ,q) . \n") ; 
To be on the safe side, repeat the test.\n"); 



■■ 3 or test [2] = 4 then 

In short, the test found an involution which behaves as an\n"); 

involution of type t_-[n-l}- in S0(2n+1 ,q) . \n") ; 

this is not a very conclusive answer; to be on the safe side,\n"); 
repeat the test.\n"); 



• 5 or test [2] = 6 then 

In short, the test found ain involution which behaves as aoiXn"); 
involution of type t_n in S0(2n+l,q) An") ; 

This is not a very conclusive answer; to be on the safe side,\n"); 
repeat the test.\n"); 



■■ 7 or test [2] = 8 or test [2] = 10 then 

The test found an involution j which behaves either as an involution\n") ; 

of type t_l in PSp(2n,q), or as an involution of type t_n in\n"); 

Omega (2n+l ,q) ; further distinction was made by not being able to find\n"); 

a 2-element from the component of type SL(2,q) which should be present\n"); 

in the centraliser C_G(j) in the sympelctic caseAn"); 

To be on the safe side, repeat the test An"); 



■■ 9 then 

The test found an involution j which behaves either as an involution\n") ; 
of type t_l = t_-[n-l3- in PSp(2n,q), or as an involution of type t_n in\n"); 
Omega ( 2n+ 1 ,q) ; further distinction was made by being able to find\n"); 
a 2-element in C_G(j) which behaves as if it belongs to\n") ; 
the component of type SL(2,q) which should be present \n"); 
in the centraliser C_G(j) in the sjnnplectic caseAn"); 
To be on the safe side, repeat the test An"); 



■■ and test [2] = 13 then 

Since n = 3, the test is much simpler: we found an involution j\n"); 
generated by a good element and such that the primitive prime divisor rank\n"); 
of some raindom product j'x * j~y of conjugates of j is at least 3. \n"); 
This DEFINITELY identifies the group as symplectic An") ; 



■■ 1 and test [2] = 13 then 

Since n = 3, the test is much simpler: we found an involution j\n"); 
generated by a good element with the primitive prime divisor ranksXn"); 

of a random products of conjugates j~x * j~y at most 2. \n") ; 
This almost definitely identifies the group as orthogonal An" ) ; 
To be on the safe side, repeat the test An"); 



if test[l] = and test [2] = 14 then 

Print(" The programme found an involution which behaves as \n"); 
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PrintC involution of type t_n in PSp(2n,q) . \n"); 

PrintC This is a DEFINITE answer given by \n"); 

PrintC a plain (without shortcuts) version of the algorithm. \n") ; 



if test[l] = 1 and test [2] = 14 then 

PrintC The algorithm has not found an involution which behaves as \n"); 
PrintC involution of type t_n in PSp(2n,q). \n"); 
PrintC This is NOT a definite answer, and it was given by \n"); 
PrintC a plain without shortcuts) version of the algorithm. \n" ) ; 
PrintC To be on the safe side, repeat the test.\n"); 

fi; 



if 1 > 4 then 

PrintC The first run of the test has not produced the conclusive result, \n"); 
PrintC the test was repeated ", test[l], " times. \n"); 

PrintC The above conclusion is generated at the last run of the test.\n"); 

fi; 

Print ("\n") ; 

PrintC In the lists Test(G) and PlainTest (G) , the first position is or l,\n"); 
PrintC according to the answer\n"); 

PrintC *symplectic* or *orthogonal* given by the test, the second position is\n") ; 
PrintC the rule number used in this decision, the third position is the actual type\n"); 
PrintC of the involution used in the identification of the group; \n"); 
PrintC (the latter was only used for the debugging of the algorithm) . \n") ; 
PrintC The fourth number is the numbers of re-runs of the test\n"); 
PrintC when the first run was inconclusive. \n") ; 

end; 



